Thermalization of magnons in yttrium-iron garnet: 
nonequilibrium functional renormalization group approach 



Johannes Hick, 1 Thomas Kloss, 2 and Peter Kopietz 1 

^^Institut fiir Theoretische Physik, Universitat Frankfurt, 
Max-von-Laue Strasse 1, 60438 Frankfurt, Germany 
2 Laboratoire de Physique et Modelisation des Milieux Condenses, 
CNRS and Universite Joseph Fourier, 25 Avenue des Martyrs, 38042 Grenoble, France 

(Dated: June 28, 2012) 

Using a nonequilibrium functional renormalization group (FRG) approach we calculate the time 
evolution of the momentum distribution of a magnon gas in contact with a thermal phonon bath. 
As a cutoff for the FRG procedure we use a hybridization parameter A giving rise to an artificial 
damping of the phonons. Within our truncation of the FRG flow equations the time evolution of the 
magnon distribution is obtained from a rate equation involving cutoff-dependent nonequilibrium self- 
energies, which in turn satisfy FRG flow equations depending on cutoff-dependent transition rates. 
Our approach goes beyond the Born collision approximation and takes the feedback of the magnons 
on the phonons into account. We use our method to calculate the thermalization of a quasi two- 
dimensional magnon gas in the magnetic insulator yttrium-iron garnet after a highly excited initial 
state has been generated by an external microwave field. We obtain good agreement with recent 
experiments. 
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I. INTRODUCTION 

Calculating the nonequilibrium time evolution of in- 
teracting many-body systems is a challenging problem 
which usually requires advanced many body methods and 
serious numerical calculations. This problem is of interest 
in many different fields^, including ultracold atoms^, col- 
lision experiments of heavy nuclei^, exciton-polariton sys- 
tems in semiconductor microcavities^, dissipative quan- 
tum syst ems ^ and pumped magnon gases in magnetic 
insulator^HI. The nonequilibrium dynamics of weakly 
interacting non-relativistic bosons has been successfully 
modeled by combining the Gross-Pitaevskii equation for 
the condensate with th e B oltzmann equation for the 
quasi-particle excitations^^. However, the collision in- 
tegral in the Boltzmann equation is usually calculated 
only perturbatively to second order in the bare inter- 
action. This approximation breaks down for strong in- 
teractions or for high densities, so that it is important 
to develop non-perturbative methods. For ./V-component 
relativistic scalar field theories Berges and co-authors^ 
have recently shown by means of a two-particle irre- 
ducible resummed large- N expansion that vertex correc- 
tions (which are neglected in the Boltzmann equation) 
qualitatively change the scaling behavior of the momen- 
tum distribution for small momenta. 

An alternative non-perturbative method to deal with 
strongly correlated systems out of equilibrium is based on 
the functional renormalization group (FRGjnma R e - 
cently, several authors have used nonequilibrium FRG 
methods to calculate the time evolution of various types 
of many body systems^^. It turns out, however, 
that the proper implementation of nonequilibrium FRG 
schemes depends on the specific problem of interest; in 
fact, many technical problems (such as the choice of the 



cutoff parameter, or the construction of sensible approx- 
imation strategies which respect causality and the con- 
servation laws) have to be solved in order to turn the 
nonequilibrium FRG into a useful and competitive cal- 
culational tool. 

In this work we shall develop a particular implementa- 
tion of the nonequilibrium FRG which is suited to study 
time-dependent nonequilibrium phenomena in Bose sys- 
tems in contact with a thermal phonon bath. To be 
specific, we identify the bosons with the magnons (i.e., 
quantized spin-waves) in a ferromagnetic insulator; how- 
ever the methods developed in this work are also useful to 
calculate the nonequilibrium time evolution of any non- 
relativistic Bose gas. Moreover, by means of a straight- 
forward modification of our method it should also be 
possible to study Fermi gases out of equilibrium. The 
Hamiltonian of our system is 



k q 



(1.1) 



where a k creates a magnon with momentum k and en- 
ergy efc, while b q creates an acoustic phonon with mo- 
mentum q and energy u> q — c\q\. The volume of the sys- 
tem is denoted by V, the operator p q = J^k a l a fc+<j re P~ 
resents Fourier components of the magnon density, and 
the magnon-phonon coupling 7 q is assumed to be propor- 
tional to y/uJq~, which is a simple consequence of the fact 
that density fluctuations associated with the creation of 
longitudinal phonons are proportional to the divergence 
of the corresponding displacement field^. 

A model Hamiltonian of the form (1.1 ) has been used 
previously by Banyai et alW^ to describe the kinetics of 
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Bose- Einstein condensation of excitons. In Sec.|TT]we will 
review the rate equation approach adopted by these au- 



thors. Here we use the model Hamiltonian (1.1 ) to study 



the nonequilibrium dynamics of magnons in thin stripes 
of the magnetic insulator yttrium-iron garnet (YIG) 21 . 
The time evolution of the momentum distribution of 
magnons in YIG has been studied experiment ally by 
means of the method of Brillouin light scatterin g ^ 6 ! 22 ! 23 !. 
Note that the effective magnon Hamiltonian of YIG is 
of the form ( 1.1 ) if we identify the operators a k with the 
creation operators of the relevant magnon band in YIGfSI. 
In the experimentally realized thin stripe geometry the 
magnons form a quasi two-dimensional weakly interact- 
ing Bose gas. If the external magnetic field H is aligned 
along the stripe axis and the angle k between the mag- 
netic field and the magnon momentum k is not too large, 
the long-wavelength energy dispersion of the experimen- 

ta lly re levant magnon band in YIG can be approximated 
bjJMUl 



(k = [h- 
x [h 



t fc 2 + A(l-/ fe )sin 2 fc ] 1 /2 



Pexfc 2 + A/fe] 



1/2 



(1.2) 



Here h — fi\H\ is the Zeeman energy associated with 
the external magnetic field (/x is the magnetic moment 
of the localized spins in YIG), p cx is the exchange spin- 
stiffness, and the energy A is the characteristic energy 
scale associated with dipole-dipole interactions in YIG. 
To describe the experiments it is sufficient to retain only 
a single magnon band and work with an effective two- 
dimensional model. The form factor f k in Eq. (1.2) can 
then be approximated by^ 



fk 



-\h\o 



\k\d 



(1.3) 



where d is the transverse thickness of the YIG stripe. As 



will be discussed in more detail in Sec. IV the parame- 
ters p ex , A and d entering the magnon dispersion (1.2) 



are all known from experiment, so that we may use our 
model Hamiltonian (1.1) to make specific predictions of 



the nonequilibrium dynamics of the magnon gas in YIG. 
Note that our Hamiltonian (1.1) does not contain di- 



rect two-body interactions between the magnons, so that 
in the absence of phonons the magnons do not interact. 
A more realistic model describing the Bose-Einstein con- 
densation of magnons in YIG should al so inc lude many- 
body interactions between the magnon a 21 * 26 l Although 
these interactions compete with the phonons to thermal- 
ize the magnon gas, we focus here on a regime where the 
dominant thermalization mechanism is due to magnon- 
phonon scattering. 



II. RATE EQUATIONS FROM THE KELDYSH 
TECHNIQUE 

The nonequilibrium dynamics of a boson Hamiltonian 



of the type (1.1) with quadratic boson dispersion e k has 



been studied in Ref. [5D] using a simple decoupling pro- 
cedure of the equations of motion, which yields the fol- 
lowing rate equation for the time-dependent momentum 
distribution n k (t) = (a k (t)a k (t)) , 

d t n k (t) = y^{[l + n k (t)}W k ,k>n h >(t) 
k' 

-[l + n h >(t)]W h ,, h n h {t)}. (2.1) 

The transition rates are given by Fermi's golden rule^ 
W k , k , = 2 7 L fe ,&(e fe - e k ,)Dl_ k ,{e k - e k .), (2.2) 
where the Bose function 

1 



b{u) 



o/3u} 



1 



(2.3) 



is the equilibrium distribution of the phonon bath at tem- 
perature T — 1/(3, and 

Dfa) = 7rsgn(w) [6(w - cu q ) + S(uj + uj q )] (2.4) 

is the spectral function of the symmetrized phonon prop- 
agator, see Eq. (C4c| in Appendix C. N ote that for 
negative frequencies the spectral function (2.4) is neg- 



ative, which is a general property of any bosonic spectral 
functional Using the fact that D„{—uS) = —Dl(cj) and 



6(-w) = -[l + 6(w)] 



b(tu) one easily verifies that 



the transition rates satisfy detailed balance 2 ^ 



W k >. k e 



(2.5) 



Numerical solutions of Eq. (2.1) for finite systems 



of bosons with quadratic dispersion can be found in 
Ref. [20]. However, for finite systems the Dirac 5- 
functions in the transition rates (2.2) lead to singular 
results. Banyai et a^PSl therefore replace the (5-functions 
(by hand) by Lorentzians of finite width, thus introducing 
an additional phenomenological parameter into the prob- 
lem. This is somewhat unsatisfactory, because the final 
results can be sensitive to the width of the broadened 5- 
functions. Here we present an alternative solution to this 
problem: within our FRG approach with hybridization 
cutoff it is straightforward to take the physical broaden- 
ing of the phonon propagators due to the coupling to the 
magnon system self-consistently into account. 

To obtain a FRG generalization of the rate equation 



(2.1), an approach based on the decoupling of the equa- 



tions of motion is not useful; instead, the nonequilibrium 
problem should be formulated in terms of self-energies 
which are obtained via functional integral techniques. 
Before developing the nonequilibrium FRG approach, 
it is therefore instructive to carefully re-derive the rate 
equation (2.1) using the functional integral formulation 
of the Keldysh methocP22!. Let us therefore recall that 
in the Keldysh method the nonequilibrium distribution 
function n k (t) is related to the Keldysh Green function 
at equal times, 



*G£(t,t) = l + 2n h (t). 



(2.6) 



3 



In order to calculate this, one has to consider G%(t,t') 
together with the retarded and the advanced Green func- 
tions, 

t G§(t 7 t / ) = e(t-t r )([a k (t),al(t / )}), (2.7a) 
iG£(t,t') = -e(t / -t)({a k (t),al(t>)}), (2.7b) 
iG«(t,t') = {{a k (t),al(t')}), (2.7c) 

where [ , ] is the commutator and { , } is the anticom- 
mutator. The corresponding phonon Green functions are 

iF?(t,t') = Q(t-t')([b q (t),bl(t')}), (2.8a) 
iFf(t,t') = -Q(t'-t)([b q (t),bl(t')]), (2.8b) 
iF q x (t ) t') = ({b q (t),b q (t')}). (2.8c) 



Since the magnons in Eq. ( 1.1 1 couple only to the hermi- 
tian combination 



fct 



(2.9) 



of the phonon fields, it is convenient to define 

iD^(t,t') = e(t-t , )([X q (t),X- q (t')]), (2.10a) 

iD*(tX) = -e(t' - t){[x q (t),x_ q {t>)]), (2.10b) 



iD«(t,t') = ({X q (t),X_ q (t')}). 



(2.10c) 



Following Kamenev^ESI, we represent the above Green 
functions as functional averages, 



iGftt.f) = {a c k {t)a k ! {t l ))^ l G^{t,t') 1 (2.11a) 

Ktf(t,f) = (a2(t)a^(t'))^ l G2 C (t 1 t'), (2.11b) 

iG%{t,t') = (a c k (t)a c k {t')) = iG c k c {tX), (2.11c) 
and similarly for the phonon Green functions, 

iif (i,f) = {b c q {t)bQ{t')) = iFCQ(t,t% (2.12a) 

%F£W) = (bQ{t)b c q {t')) = iF« c (*,*'), (2.12b) 

tff (f,f) = (b c q {t)b c q {t')) = iF^ c {t,t'). (2.12c) 



The functional averages are defined as 



4(t)a X k (t>)) = / V[a,aXby s ^ a ^al{t)a x k (t'), 



(2-13) 

and similarly for (b q (t)b q (t')). The superscripts A, A' e 
{C, Q} label the classical (C) and quantum (Q) com- 
ponents of the fields, which are related to the contour 
components a k (t) and b k (t) via 



<**(*) = 


i 

71 


[<£(*) 


+ «£(*)] , 


(2.14a) 


<£(*) = 


i 

71 


KW 


-<£(*)]> 


(2.14b) 


= 


i 

71 


[&£(*)- 


f 6fc(*)], 


(2.14c) 


#(*) = 


i 

71 


[&£(*)- 


-b k {t)}. 


(2.14d) 



to 



FIG. 1. (Color online) The Keldysh contour connects some 
initial time to slightly above the real axis with some final time 
t > which is larger than all times of interest. It then crosses 
the real axis and goes back to to slightly below the real axis. 
In accordance with the definition of time-ordering we draw 
later times to the left. 



Here the superscripts ± attached to the fields a k (t) and 
&k (t) refer to the upper and lower branch of the Keldysh 
contour shown in Fig. [I] The Keldysh action is given by 



i, b, b] = 



dt J dt'l 5>£(*),a2(t)) 

+£(&£(*),&?(*)) 

q 



o (c^r 1 

(i^)" 1 
(F «)-i -(F*)-iF K (Ff)-i 



Q 

fe+q 




(2.15) 



Here we have defined the phonon fields 

X x (t)^^.[b x q (t) + b x _ q (t)] , X = C,Q, (2.16) 



and Gq,Gq, and Gq as well as F^F^, and F^ are 
infinite matrices in the time-labels (we omit for simplic- 
ity the momentum label) , with matrix elements given by 
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< < 



a c a c 



a c (ft 



< < • < < — < 



:x c <x c ex® <x<* 

FIG. 2. (Color online) Graphical representatio n of t he bare 
magnon-phonon vertices of the Keldysh action (2.151. Solid 
arrows pointing into the vertices represent aft , while outgo- 
ing solid arrows represent aft . Dotted arrows represent the 
corresponding quantum components aft and aft , Wavy solid 
lines represent the classical component X c of the phonon field 
defined in Eq. (2.161, while wavy dotted lines represent the 



quantum component X 



where go and fo are diagonal matrices in the time labels 
which contain the magnon- and the phonon distribution 
functions in the non- interacting limit, 



Mh,w = S{t - t')g , k = S(t - t')[l + 2(a{a k ) }, (2.18a) 
[/o]„tf = S(t - t')/o,g = 6(t - t')[l + 2{b\b q ) }. (2.18b) 



the Green functions (2.11a 2.12c) in the limit where the 
magnon-phonon coupling is switched off. A graphical 
representation of the interaction vertices in Eq. ( |2.15 1 is 
shown in Fig. [2] Note that the lower diagonal blocks of 
the inverse Gaussian propagators in Eq. (2.15) are actu- 



ally infinitesimal regularizations of the continuous func- 
tional integral, 



K 



o (Go 
(F A 



A\-l 



^0 



= 2«?7/o, 



(2.17a) 
(2.17b) 



Note that we restrict ourselves to initial states with- 
out correlations. Therefore the initial density matrix 
appears only via the initial conditions for the indepen- 
dent one- and two-point Green's function^HHHH Since 
the Keldysh action (2.15) is a Gaussian functional of the 



phonon fields, we may perform the integration over the 
phonon fields exactly and obtain the following effective 
Keldysh action for the magnon fields, 



S[a,a} = Jdtj dt'J2(a%(t),a%(t)) 



(G6 4 ) 



A\-l 



(G*)" 1 -(G R )-^(G A )^ ) tv \*W) 



Jdtj di'i^ 7 ^< q (i,<')p C ?W[p^(i') + ^ Q (^)] +D A Jt,t')[p™(t) + p^(t)]pQ c (t') 



(2.19) 



where DQ q (t,t'), X = R,A,K, denote the non- 
interacting limit of the symmetrized phonon propagators 
defined in Eqs. ( |2.10a -2.10c), and we have defined 



Pf 



(*) = X>fc(*)0M.,(t), A,A'e{G,Q}. (2.20) 



It is convenient to collect the retarded (G R ), advanced 
(G A ) and Keldysh (G K ) components of the Green func- 
tion into a matrix Green function, 



G = 



\G] CC \G] CQ 







G K G R 



G- 







whose inverse has the block structure 

c -i = ( {G A )~ l 

1 (G 11 )- 1 -(G B )- 1 G K (G A )- 1 



(2.21) 



(2.22) 



Defining the nonequilibrium self-energy via the matrix 
Dyson equation, 



G- 1 = Gg" 1 



(2.23) 



we have in the Keldysh basis, 



£ = 



[S] CQ 
[S1 QC [£1 QQ 





( 




)= 


{ E R 


t K ) 



with 



.[ G -i]QQ = {G R )- 1 G K {G A Y 



(2.24) 



(2.25) 



By taking appropriate matrix elements of the nonequi- 
librium Dyson equation (2.23), we obtain a quantum ki- 
netic equation for the distribution function. Actually, 
as reviewed in Appendix A, one can rewrite the kinetic 
equation in several different forms, depending on the 
choice of basis and on the re-shuffling of the terms in 
the Dyson equations. For our purpose, it is most conve- 
nient to parametrize the Keldysh block in terms of the 
distribution matrix g by setting 



G K = G R g^ - gG J 



(2.26) 



where in the non-interacting limit the matrix g reduces 
to the non-interacting distribution matrix go given in 
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Eq. (2. 18a |. It is then easy to show (see Appendix A) 
that the Dyson equation ( |2.23 ) implies 

-i{Mg-$M) = E in -E out , (2.27) 

where the infinite matrices M, S m and E out can be writ- 
ten as 



M = M o - £ 
E in = it K 



(2.28) 
(2.29) 

(2.30) 



Here the matrix elements of Mq in the momentum-time 
basis are 

[Mo]kt,w = S htk > [iS'(t - t') - e k 5(t - t')} , (2.31) 

where 5'(t) — 4;S(t) is the derivative of the Dirac 5- 
function. Moreover, we have introduced the mean (de- 
noted by a superscript M) and the spectral (denoted by 
a superscript I for imaginary part) combinations of the 
retarded and advanced self-energies^, 



t M = ~[t R + ± A ] , t T =i[t R 



(2.32) 



Our notation for E m and S out anticipates that these 
quantities are related to the "in-scattering" and "out- 
scattering" contributions to the collision integral of the 
kinetic equation. 

To describe the approach towards equilibrium, it is use- 
ful to express our kinetic equation in terms of Wigner 
transforms. Given any infinite matrix A in the time la- 
bels with matrix elements we define its Wigner 
transform A(t; uj) by performing a partial Fourier trans- 
formation with respect to the time difference s = t — if , 



A(r;u) 



dse luJS [A] TJl 



(2.33) 



which is a function of the average time r = and the 

The 



frequency w associated with the difference t — t'. 
Wigner transform of our kinetic equation (2.27) is 



d T Rcg k (r] w) + 2(u - e k )lmg k (T;uj) 



Vg- 



M - 



n 



M\ 



<(r,u) 



(2.34) 



of g* is given by the function g*(r;ui), and (AB)t T>u ) de- 
notes the Wigner transform of the product of two matri- 
ces in the time labels. Although in general such a prod- 
uct does not factorize into the product of Wigner trans- 
forms of A and B, there is an approximate factorization 
to leading order in a gradient expansion, see Eq. (A26|. 



Assuming that such an expansion is justified, we obtain 
to leading order 
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2£f (r; W )Im 5fc (r; W ), (2.35) 



(t T g + tft 1 ) « e£(t; W )Re fffe (r; u), 



(2.36) 



so that our kinetic equation (2.34) becomes 



d T Reg k (T:uj) + 2[cj -e k - (r; w)]Img fe (r; uj) 

= i^{r;u>) - E^r; W )Re 5fc (r; W ). (2.37) 

Next, let us neglect the term involving the imaginary part 
of g k (r\uj) and set uj = e k . Defining the quasi-particle 
distribution function g k {r) = Reg k (T;uj = e k ) we arrive 
at the simplified kinetic equation 

d T g k (r) = i£*(r; e k ) - ££(t; e k )g k (r). (2.38) 

Setting g k (r) = 1 + 2n k (r), this equation can also be 
written as 



d T 2n k (r) = Ejftr) 



^out 



(r), 



(2.39) 



where the "in-scattering" and "out-scattering" self- 
energies are given by 

4 n (T)=zSf(T; efc ), (2.40a) 
S^ t (r) = S^r;e fe )[l + 2n fe (r)]. (2.40b) 

To make further progress, we need approximate expres- 
sions for the nonequilibrium self-energies Yj^(r\€ k ) and 
E^(r;efc). In Appendix B we give perturbative ex- 
pressions for these self-energies to second order in the 
magnon-phonon coupling 7 q . Note that the simplest way 
to obtain the self-energies to this order is to perform a 
Hartree-Fock decoupling of the interaction term of the 
effective boson action (2.19). Within the same approx- 



imation as above (factorization of Wigner transforms, 
neglecting renormalization of the magnon energies, see 
Appendix C for details) we obtain 

iEffae*) = ^-^2{[l + n k/ (T)]W k ^ k + W k ^n k ,(T)}, 



(2.41) 

Sfc(r;e fc ) = i ^{[1 + n k . {r)]W kl , k - W k , k ,n k ,{ T )}, 
fc' 

(2.42) 



where the transition rates W kyk ' are given in Eq. (2.2) 



Subst ituting Eqs. (2.41 2.42) into the kinetic equation 
(2.38) and renaming r->(we finally arrive at the rate 



where we have used the fact that the Wigner transform equation (|2T| 



For finite systems the (5-functions appearing in the 
phonon spectral function (2.4) appearing in the pertur- 



bative expression (2.2) for the transition rates should be 
regularized. The simplest possibility is to perform an 
ad hoc replacement of the (^-functions by Lorentzians 
with finite width 20 , but this procedure introduces an 
additional phcnomcnological parameter into the prob- 
lem. Microscopically, the coupling of the phonons to the 
magnon system automatically generates such a broaden- 
ing. We therefore introduce renormalized phonon propa- 
gators involving the phonon self-energies (t;uj), where 
X = R,A, K. Defining the mean and the spectral 



G 



<....{ n«i 



' [U A h-< 



G R G A ' 



^Or< + .3-* 

qR qA 



FIG. 3. (Color online) Diagrams representing the phonon 
self-energies Tl H , Tl A and Ii K to second order in the magnon- 
phonon interaction. Wavy lines represent the external phonon 
fields b c and b Q , while arrows represent the magnon Green 
functions G R , G A and G K . We use the convention that the 
classical field components are represented by solid lines, while 
the quantum components are represented by dotted lines, see 
also the caption of Fig. [2] If we introduce in our FRG ap- 
proach a hybridization cutoff only in the phonon propagators, 
these diagrams are finite at the initial cutoff scale. 



combinations of the retarded and advanced phonon self- 
energies, 

jjM _ ijjjfl + jp4] ; fii _ i^jR _ jjA] ( 2 43 ) 

the Wigner transforms of the renormalized phonon prop- 
agators can be written as 



-II^(t;cj)' 
1 

u-u q - n^(r;w)' 



(2.44a) 
(2.44b) 



n^(r; W )]^+[inj(r; W )]' 



i F q(i",u) = coth F^(t;u) 



j3u 



(2.44c) 
(2.44d) 



The diagrams of the phonon self-energies up to second 
order in the magnon-phonon coupling "/ q are shown in 
Fig. [3| Approximating the Wigner transforms as dis- 
cussed above, and neglecting the real parts of the self- 
energies, we have H M s» and II ^ ps ^H 1 rj — T1 A , where 
the spectral component of the phonon self-energy is 



2tt 



U l ( T ' W ) = Y ^1 6 + 6k -1 ~ 6fe ) 
k 

x {n k _ q (t) - n k (t)} . 



To obtain regular results also for finite magnon gases, we 
evaluate the right-hand side of Eq. (2.45) in the limit of 
infinite volume, where IT^ (r; u) is a continuous function 
of all of its arguments. Let us emphasize that we take 
the limit V — > oo only in the phonon self-energies, so 
that with the resulting regular transition rates we can 
still use our rate equation (2.1) to describe the kinetics 



of a finite magnon system with discrete spectrum. 



III. FUNCTIONAL RENORMALIZATION 
GROUP APPROACH 

Formally exact FRG flow equations for the irre- 
ducible ver tices of any bosonic many-body systems out of 
equilibriu m 1 16 * 17 ! can be obtained as a special case of the 
general hierarchy of FRG flow equations^ which follows 
from the Wetterich equation^. While it is straightfor- 
ward to write down these equations, the proper identifica- 
tion of a sensible cutoff procedure which does not violate 
causality and respects the fluctuation-dissipation theo- 
rem under equilibrium conditions is rather delicate prob- 
lem. It turns out that for the coupled magnon-phonon 
system considered here a modified hybridization cutoff 
scheme^^ is most convenient. In this scheme the in- 
finitesimal imaginary parts ±«ry defining the boundary 
conditions of the retarded and advanced propagators are 
replaced by finite quantities ±zA, where A is identified 
with the flow-parameter of the renormalization group 
procedure. Moreover, in the hybridization cutoff scheme 
the same replacement is also made in the Keldysh blocks 



of the inverse Gaussian propagators, see Eqs. (2.17a) and 
( 2.17b[ ). For bosonic systems, however, one shou d be 
careful to take into account that the spectral function 
should be positive for positive frequencies, and negative 
for negative frequencies^. The imaginary part of the 
spectral component of the self-energy therefore changes 
sign at zero frequency. Hence, in the bosonic version of 
the hybridization cutoff scheme the infinitesimal imag- 
inary part r\ should be replaced by r\ — > Asgn(w). In 
contrast, for fermions the spectral function is positive for 
all frequencies, so that the hybridization cutoff should bc 
introduced by replacing r\ — \ A. 

In this work, we introduce a hybridization cutoff only 
in the phonon propagators, while keeping the ±ir] in the 
magnon propagators infinitesimal. In this way we do 
not artificially smear out the magnon distribution, so 
that we can still describe the nonequilibrium dynamics 
in terms of a kinetic equation describing the evolution 
of the occupation probabilities of sharp energy levels. 
The cutoff-dependent Keldysh action S\[a,a,b,b] with 
a hybridization cutoff in the phonon propagators c an be 
obtained from the action S[a,a, b,b] given in Eq. (2.15) 



(2.45) 



by replacing the Gaussian phonon propagators by cutoff- 
dependent propagators, 
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S A [a,a,b,b] = dt dt'l £>£(t),a£(*)) 



(G^)- 1 
(G^)- 1 2ir,g 



(*)) 



TO 



<(0 



C Q 

k+q u k 



A 

,C 

i fc+q d fc 



k+q a k 



(3.1) 



Because for A -> oo all phonon propagators vanish, the 
FRG flow should be integrated with the boundary condi- 
tion that all irreducible pure magnon vertices vanish for 
A — > oo. This implies that all closed loops involving only 
magnon propagators do not vanish in this limit, leading 
to a non-trivial initial condition for the FRG flow. This 
is ana logou s to the bosonic momentum transfer cutoff 
scheme^^ in coupled fermion-boson systems, where at 
the initial scale all pure boson vertices are finite. In par- 
ticular, in the present problem the phonon self-energy 
diagrams shown in Fig. [3] are finite at the initial scale. 
Neglecting the renormalization of the real part of the 
phonon self-energies, the scale-dependent phonon propa- 
gators are in our cutoff scheme given by 



F* q {r;u>) 



iFL{r,u) 



ui - Lu q + iAsgnw + _(r;o>) 



1 



ui - Lu q - iAsgnw - |n^ 9 (r; ui) 
{u - u q y + [Asgnw + |nj (t; lj)]' 



, (3.2a) 
, (3.2b) 



(3uj 



coth[ ^ ] Fl q {T-uj), 



(3.2c) 
(3.2d) 



where the imaginary part of the scale-dependent phonon 
self-energy is approximated by an expression similar to 
Eq. |2~4o), 



2?r 



k 

x {n^k-q (t) - iiA,fe (r)} , (3.3) 



where n\ k( T ) is the scale-dependent magnon distribu- 
tion function, and it is understood that for our numerical 
calculations we replace V~ 1 '£ lk I d 2 k/(2n) 2 to ob- 
tain a smooth phonon self-energy in our two-dimensional 
system. Note that the scale parameter A can be con- 
sidered as an additional contribution to the imaginary 
part of the phonon self-energy, which arises from the hy- 
bridization between phonons and magnons. Introducing 
the hybridization function 



1, 



<-[t, r 



m 

<-{sy-< 



D R D K 
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FIG. 4. (Color online) Graphical representation of the trun- 
cated FRG flow equations for the nonequilibrium magnon self- 
energies with hybridization cutoff in the phonon propagators. 
The dots over the self-energy symbols on the left-hand side 
represent the derivative with respect to the cutoff parameter 
A. Slashed wavy lines on the right-hand side represent single- 
scale phonon propagators, wh ich can be obtained by taking 
the partial derivative of Eqs. (3.5a ]3T5d | with respect to the 
explicit A-dependence. Note that the diagrams on the right- 
hand side can be obtained from the self-energy diagrams given 
in Appendix B (see Fig. [6| by simply replacing the phonon- 
propagators by the corresponding single-scale propagators. 



the symmetrized phonon propagators can be written as 
cjq - iAq(r;cj) 



DlqH 



UJ 2 


- [cjq - iAq(r; 


uj)} 2 




Wq + iAq(r; w 


) 




- [cjq + zAq(r; 


c)] 2 




Aq(r;w) 






-Cq)2+Aq(T 


u>) 




Aq(r;w) 




( W 


+ ^q) 2 + A2(r 





iD^ q (u>) = coth 



(3uj 
~2~ 



D{,q(u)- 



(3.5a) 
(3.5b) 

(3.5c) 
(3.5d) 



Aq(r;w) = Asgnw + -n^ g (r;w), 



(3.4) 



The general form of nonequilibrium FRG flow equa- 
tions for interacting bosons in the Keldysh basis has been 
written down in Refs. [161 and I17j . Here we use a sim- 
ple truncation where the renormalization of the magnon- 
phonon vertex is neglected. The relevant FRG flow equa- 
tions for the nonequilibrium magnon self-energies with 
hybridization cutoff in the phonon propagators are shown 
graphically in Fig. [4] For the explicit evaluation of these 
expressions, we use the same approximations as in the 



derivation of the perturbative rate equation (2.1| out- 
lined in Sec. [TTJ namely: 

1. Keep only the leading terms in the gradient ex- 
pansion so that Wigner transforms of products of 
infinite matrices in the time labels factorize. 

2. Neglect the contributions in the kinetic equation 
describing the renormalization of the magnon ener- 
gies. 

3. Take the frequency argument of the self-energies on 
resonance, w — e k . 

After straightforward manipulations analogous to those 
described in Sec. [TT| we obtain a system of three cou- 
pled partial differential equations for the scale-dependent 
magnon distribution function n,A,k{ T ) an d the scale- 
dependent nonequilibrium self-energies H A k (r,e k ) and 

fc( T ' £ k)- The FRG version of the kinetic equation is 
formally identical with the perturbative kinetic equation 
(p8l 



3 T 2n A ,fc(r) = iSf \ k (r;e k 



Sl, fe (r;e fe )[l + 2n A)fe (r)]. 

(3.6) 



In contrast to the perturbative equation (2.38) where 



the nonequilibrium self-energies can be expressed 
as integrals depending on the distribution function 
[see Eqs. (|2.41| |2.42"|)], in our FRG approach the self- 



energies are determined by flow equations relating their 
derivatives with respect to the cutoff A to the A- 
dependent distribution function, 



9 A iSf ife (r;e fe ) = i + n A , k , (r)]W k , , 



y^ 1 

k' 



5aS{ fc (r;e fc ) 



+Wfc,fc'nA,fc'(r)}, 
= ^2{[l + n A ,k'(T)]W k >,k 



V 



-Wfc,fc'nA,fc'(r)}. 



(3.7) 



(3.8) 



Here W k ,k' is the partial derivative of the cutoff- 
dependent transition rate with respect to the explicit A- 
dependence, 



W k ,k> = 27fe_fe'i>A,fe-fe'(efe - e k ,)b{e k - e k ,), (3.9) 

same structure as Eq. (! 
phonon spectral function D Aq (uj) repl. 

A,(r; W )( W _ U Q 



which has the same structure as Eq. ( 2.2[ ), but with the 

aced by 



^A, q M = Trsgnw^ 



L' 



where 



L' A (u) = d A 



A 



A 2 



1 u 2 



A,(t;w) 



A 2 



it (u 2 + A 2 ) 2 



(3.10) 



(3.11) 



is the derivative of a normalized Lorentzian with respect 
to its width. Note that the right-hand side of Eq. (3.10) 



implicitly depends on the time t via the r-dependence of 



the phonon self-energy iT A q (r; w) given in Eq. ( 2.45 1 . To 
simplify the numerical evaluation of the above system of 
integro-differential equations, we shall replace in Sec. |TV| 
the phonon self-energy by its equilibrium value at van- 
ishing cutoff, which is approached for t — > oo. In other 
words, we replace 



(3.12) 



This procedure can be justified by noting that for suffi- 
ciently large values of the cutoff parameter A the contri- 
bution from the phonon self-energy to the hybridization 
function Aq(r;w) defined in Eq. ( |3.4[ ) can be neglected, 
so that the only for small A the effect of the phonon 
self-energy becomes important; but in this regime the 
magnon distribution is already close to its equilibrium 
value. 

Unfortunately, the above system of equations violates 
particle number conservation, in contrast to the pertur- 
bative kinetic equation (12. II), which obviously implies 



d_ 

dt 



J2n k (t)=0. 



(3.13) 



The fact that Eqs. (3.6 3.8) violate particle number con- 



servation is related to our approximation of neglecting 
the renormalization of the real part of the magnon ener- 
gies. Note that particle number conservation is related to 
the U (l )-sy mmctry of the bare action S A [a, a,b,b] given 
in Eq. (3.1). To see this more clearly, let us multiply 
the magnon fields in Eq. ( |3.1[ ) by an arbitrary scale- and 
time-dependent phase factor (X = C,Q), 



«fc (*) 



- J°<A(t) n fX 



- p -ia A .{t) Ji tX 



«fe (*)■ (3-14) 



Then the Keldysh action (3.1) changes as follows, 

S A [e~ ia a', e ia a', b, b] = S A [a', a', b, b] 
da A {t) y-^ r , c 



dt- 



dt 



k 



[d' k a (t)a^(t) + a' k Q (t)a / k c (t)].(3.l5) 



The additional term can be absorbed via a shift in the 
retarded and advanced self-energies, 



R \«+5(t-t>) daA{t) 



dt 



(3.16a) 
(3.16b) 



The extra terms drop out in the spectral self- 
energy T, 1 = — fl A ), but survive in the average 
S M = ^[t R + By completely neglecting the sym- 

metric combination S M in our derivation of Eq. (3.6), 



we have lost the possibility of explicitly implementing the 
condition of particle number conservation. To impose the 
constraint of constant particle number on the FRG flow, 
we add the corresponding symmetric self-energy E M to 
the left-hand side of the kinetic equation. For our pur- 
pose, it is sufficient to approximate 



-S(t-t')fi A (t). 



(3.17) 
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The only effect of /i A (i) on the kinetic equation (3.6) is 



to replace the derivative d T by the corresponding covari- 
ant derivative d T + ^a(t). Note that the time-dependent 
Lagrange multiplier /xa(t) enforcing the constraint of 
particle-number conservation appears in the quan tum dy- 
namics as a gauge -fielcP. As a result, Eq. (Eol) should 
be replaced by 

[d T + ^ A (T)]2n A>fe (r) 
= i^ k (r; e k ) - S Ajfc (r; e fe ) [1 + 2n A , fe (r)] . (3.18) 

The energy ^a(t) plays the role of a scale- and time- 
dependent Lagrange multiplier, which should be adjusted 
such that the total particle number TV = X)fc n A. is 
conserved for all times r and for all values of the cutoff 
A, implying 



(3.19) 



Substituting Eq. (3.18) into the right-hand side of this 



equation, we conclude that the scale- and time dependent 
Lagrange multiplier is explicitly given by 

Ma(t) = ^^{i^'lki^ek) 

k 

-Ei )fc (r;e fc )[l + 2n A , fc (r)]}. (3.20) 



Combining Eqs. (3.18) and (3.20) with Eqs. (3.7), (3.8) 



and (2.45) we obtain a closed system of particle-number 



conserving FRG flow equations which can be used to cal- 
culate the nonequilibrium time evolution of the magnon 
distribution function. In the following section we shall 
use these equations to calculate the thermalization of the 
magnon gas in YIG. 



IV. THERMALIZATION OF MAGNONS IN YIG 

Quantized spin-waves (magnons) in ordered magnets 
obey to a very good approximation Bose statistics, so 
that magnetic insulators can serve as relatively easily ac- 
cessible model systems for investigating the nonequilib- 
rium dynamics of bosons. A particularly well character- 
ized magnet is yttrium-iron garnet (YIG) 21 . The mo- 
mentum distribution function of the magnon gas in YIG 
has been probed exp eri ment ally using the technique of 
Brillouin light scattering^ 22 ' 23 ' 3 ^ . 

At long wavelengths the energy dispersion of the ex- 
perimentally relevant magnon band in thin YIG-stripes 
can be approximated by Eq. (1.2). The values of the 
exchange stiff ness p cx and the dipolar energy scale A 
entering Eq. (1.2) can be written as p cx = JSa 2 and 



A = 47r/iA/ s , where J « 1.29 K is the nearest neighbor 
exchange constant of YIG and the saturation magneti- 
zation is given by AttM s ps 1750 G. Here a 12.376 A 
is the lattice constant and \x = g/iB is the magnetic mo- 
ment of the localized spins^. If we arbitrarily set the 



effective g-factor equal to two so that fi — 2fig, then 
the effective spin is S = M s a 3 //1 w 14.2. Due to an in- 
terplay between exchange interactions, dipole-dipole in- 
teractions, and finite-size effects, the magnon dispersion 
of YIG exhibits two generate minima at finite momenta 
±A: m i n £, where z is the direction of the external mag- 
netic field. Having fixed the energy dispersion, the only 
adjustable parameters of our model Hamiltonian ( 1.1 ) are 
the phonon velocity c = w q /|q| and the magnon-phonon 
coupling constant 7 = y q / \/\q\- 

Highly excited nonequilibrium states of the magnon 
gas in YIG can be generated via parametric pumping 
with a time-dependent external microwave source. Once 
the microwave source has been switched off, the thermal- 
ization of the magnon gas can be directly observe d 22 * 23 !. 
In particular, in Ref. |23j the time evolution of the mo- 
mentum distribution function (t) of the magnon gas in 
YIG has been measured. The data presented by Demi- 
dov et aZP^show how a highly excited initial state evolves 
towards a thermal equilibrium state where the occupa- 
tion of the magnon modes at the bottom of the spec- 
trum is strongly enhanced. Whether or not such a state 
should be called a Bose-Einstein conde nsate of magnons 
seems to be a m atter of semantic d 39 * 40 l We have argued 
previousl} * 26 * 41 * that the experimentally observed strong 
enhancement of the magnon distribution is not accompa- 
nied by superfluidity, because a macroscopic occupation 
of a certain magnon mode is simply equivalent with a 
change in magnetic ordeJ^. 

In Fig. [5] we show our results for the time evolu- 
tion of the magnon distribution obtained from the nu- 



merical solution of the FRG rate equations (3.7 3.8) 



and the particle number conserving kinetic equation 
given in Eqs. (3.18 3.20). For the numerical calcula- 



tions we have used the energy dispersion (1.2) with ex- 



perimentally relevant parameters discussed above and 
stripe thickness d = 5 /im. The minima of the dis- 
persion are located at \kn\ = 5.5 x 10 4 cm _1 = fc min and 
k± = 0. For simplicity we consider a square stripe with 
the length of L = 4.58 /im and work with a trun- 
cated momentum space such that |fcii| < 41.1 x 10 4 cm" 1 
and |fcjj < 13.7 x 10 4 cm _1 , corresponding to 61 x 21 
points in momentum space. The rate equation and the 
flow equations are solved simultaneously as described in 
Ref. [T7]. In order to calculate the imaginary part of 
the phonon self-energy as discussed in Sees. |H] and |III| 
we have used a two-dimensional adaption of the tetrahe- 
dron method^. The initial cutoff Ao was chosen to be 
404 J7q, where r2 = ck mm is the relevant energy scale for 
the phonon dynamics. We have set the phonon velocity 
in YIG to c = 7.15 x 10 5 cm/s as measured in Ref. |44j . 
For the numerical calculations we have divided the time 
interval into 3240 equidistant steps up to a final time of 
670 ns. For small values of the flow parameter A < ilo 
the dependence of the magon self-energies on A is rather 
strong. We have therefore divided the range of A into two 
intervals, [0,4r2 ] and [4f2 ,Ao], each equidistantly dis- 
cretized with 500 points. Our theoretical results shown 
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FIG. 5. (Color online) These contour plots show an inter- 
polation of the magnon distribution of YIG on a 7 x 7-mesh 
in momentum space for different times. The white contour 
lines denote the magnon spectrum of YIG for fixed energies. 
The crosses mark the minimum of the dispersion. The re- 
sults are presented such that they are directly comparable 
with the experimental data of Ref. [23] . We have shifted the 
zero of time to the moment where the microwave pumping is 
switched off. For the calculation of the dispersion of YIG we 
have used the parameters from Ref. |24| . which differ slightly 
from the parameters used in Ref. [23] . The initial state is cho- 
sen such that magnons with an energy of 4.1 GHz in the area 
|&x| < 0.5fc m i n are equally occupied, which resembles the ini- 
tial magnon distribution in the experiment after the pumping 
has been switched off. The strength of the magnon-phonon 
interaction has been adjusted such that the occupation of the 
ground state is saturated after approximately 670 ns, which 
gives 7 = 2.0 x 10~ 7 cmVGTTz". Finally, we have used typical 
experimental values for temperature (T = 300 K), magnetic 
field (H = 1000 Oe), and magnon density (4.77 x 10 6 /^m 2 ). 



in Fig. [5] agree quite well with the corresponding exper- 
imental data shown in Fig. 3 of Ref. 23J. Note that in 
our simple model the order parameter (a^(t)) vanishes, 
so that the global U(l) symmetry of the Hamiltonian 
( 1.1 ) is not spontaneously broken and our magnon gas is 



not superfluid. Nevertheless, our model describes the ex- 
perimentally observed thermalization towards a magnon 
state with strongly enhanced occupation at the bottom 
of the spectrum, indicating that the experiments observ- 
ing BEC of magnons in YIG can be explained without 
postulating that the magnon state is superfluid. 



Of course, the thermalization of the magnon gas in 
YIG is not exclusively driven by magnon-phonon inter- 
actions, but also by many-body scattering processes in- 
volving two or more magnons. In particular, interactions 
involving three and four powers of the magnon opera- 
tors should be added to obtain a more realistic model 
for the magnon gas in YIGP. While the three -magnon 
vertices yield the dominant contribution to the magnon 
damping^, the four-magnon vertices compete with the 
magnon-phonon interaction to re-distribute the excited 
initial state over the magnon spectrum. The inclusion of 
these scattering processes into our FRG formalism is in 
principle possible but requires substantial extensions of 
the present formalism which are beyond the scope of this 
work. 

Finally, let us point out that the magnon-phonon inter- 
action in our model Hamiltonian ( 1.1 1 conserves the num- 



ber of magnons because the phonon field X q couples only 
to the magnon density p q . However, we expect that in 
YIG there should exist another type of magnon-phonon 
interaction involving the combinations a k a)_ k _ q X q and 
a-kCLk+ q X_ q . This type of magnon-phonon coupling 
does not conserve the magnon number and leads to an ef- 
ficient energy transfer from the magnons into the phonon 
system. Microscopically, this magnon-phonon interaction 
arises from the fact that the quasi-particle operators a k 
describing the physical magnon excitations of YIG are re- 
lated via a canonical (Bogoliubov) transformation to the 
Holstein-Primakoff bosons resulting from the bosoniza- 
tion of the spin operators; see Ref. |46j for a microscopic 
derivation of these terms for a frustrated antifcrromag- 
net. Even if the couplings of these processes in YIG are 
small, we expect that these processes are important to 
describe the experimentally observecP long-time decay 
of the BEC. 



V. CONCLUSIONS 

In summary, we have studied the thermalization of the 
magnon gas due to coupling to the thermal phonon bath 
in thin films of the magnetic insulator yttrium-iron gar- 
net. Although we have only retained the interaction pro- 
cesses which conserve the total number of magnons, our 
theoretical predictions for the time evolution of the mo- 
mentum distribution of the magnons agrees quite well 
with experimental results. In particular, the accumula- 
tion of magnons at the bottom of the spectrum is cor- 
rectly described within our approach. 

We have intentionally considered in this work only the 
case where all nonequilibrium averages (a k (t)) vanish 
identically^. These averages play the role of order pa- 
rameters; keeping in mind that the magnon operators 
in YIG describe the quantized fluctuations around the 
classical ground states, macroscopic values of (a k (t)} for 
one or several wave-vectors describe a macroscopic re- 
organization of the magnetic stated To describe this 
phenomenon microscopically, it is necessary to take the 
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magnon-magnon interactions into account. 

The present work describes also some technical ad- 
vance: we have developed an implementation of the 
nonequilibrium functional renormalization group method 
which clearly exhibits the connection to a conventional 
rate equation approach. In a sense, our FRG approach 
can be considered as a simple renormalization group ex- 
tension of a rate equation approach. The main advan- 
tages of our approach are that our FRG resummation 
goes beyond the simple second order approximation for 
the collision terms in the kinetic equation, and that our 
method allows us to include in a simple way the feed- 
back of the magnon dynamics on the phonon propagators 
which in turn determine the transition rates between the 
magnon states in the rate equation. We have also devel- 
oped a generalized hybridization cutoff scheme for inter- 
acting bosons. It is possible to generalize our approach to 
include magnon-magnon interactions or magnon number 
violating magnon-phonon interactions. 
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APPENDIX A: QUANTUM KINETIC 
EQUATIONS 

In this appendix we review several possibilities of writ- 
ing down quantum kinetic equations. Although in the 
main text we only make use of the Wigner transformed 
kinetic equation for the distribution function gk(T] ui) de- 
fined via Eq. (2.26), it is instructive to work out the 



precise relation between our parametrization and a for- 
mulation of the nonequilibrium problem in terms of the 
two-time Keldysh Green function. 



1. Two-time Keldysh Green function 



Writing the nonequilibrium Dyson equation (2.23) as 
(G^ 1 - S) G = I (Al) 
we obtain following three equations for the sub-blocks, 
{G R )- X G R = 1 + £ R G R , (A2a) 



(Gq)~~ 1 g a = i + Tj A G A . 

(G R )- 1 G K = £ K G A + £ R G K , 



(A2b) 
(A2c) 



where I is the unit matrix in the momentum- and time 
labels. Alternatively we can also consider the correspond- 
ing "right Dyson equation" , 



G (Gg 1 - S) = I, 
which implies the following relations, 

d R (G R y 1 = i + 6 R £ R , 



G (Gq)~ = I + G £ 
G K (G A y 1 = G R £ K + G K t A . 



(A3) 

(A4a) 
(A4b) 
(A4c) 



Subtracting the Keldysh component of the left and right- 
hand sides of the Dyson equations ( A2c|A4c l, we obtain 
the kinetic equation 

[M , G K ] = t R G K - G R £ K + t K G A - G K ± A , (A5) 

where [ , ] denotes the commutator, and M is defined 
in Eq. <|2.31b. Note that 



(G R )- l = M + ir,I , (Go 4 )- 1 = M Q - i v I, (A6) 



A\-l 



and 



implying 



\MqG k ]u' = (id t -e k )G K (t,t'), 



\tt> = you t - tk)<^ yu, i ), (A7a) 
[G K M ]w = HcV - e k )G K (t,t'), (A7b) 



[M ,G K ]w = {id t + id t ,)G K {t,t') 



(A8) 



We conclude that in the time domain the quantum kinetic 
equation (|A5|) assumes the form 



(id t + id v )G K {t,t') 



dt^viM^ihX) - G R (t,h)i: K (h,f)} 



+ / dt 1 [i: K (t,t 1 )G A (t 1 ,t')-G K (t,t 1 )i: A (t 1 ,t')] 



(A9) 

To identify the operators which reduce in the appropriate 
limit to the collision integral in the Boltzmann equation, 
it is convenient to introduce^! £ M = ^[t, R + t, A ] and 



S J = i[T, H - see Eq. (|2_32j). The inverse relations 
are 



t M - -t 1 
2 



t A = t M + -t 1 . 



(A10) 



A similar decomposition is also introduced for the re- 
tarded and advanced Green functions, 

G M = ^[G R + G A ] , G 1 = i[G R — G A ], (All) 



so that 



G R = G m_ 1 _qI G A = G M + -G I . (A12) 
2 2 K ' 



12 



Defining M = M - E M [see Eq. (2.28)], we may write 



the kinetic equation ( A5 ) as 



M,G K ] - [E^,G M ] =C 



C° 



where we have introduced the collision matrices 



G in = ^{E K , G 1 }, 



(A13) 

(A14) 
(A15) 



Recall that { , } denotes the anticommutator. The no- 
tation indicates that after the usual approximations the 
matrix elements of of G ln and G out can be identified with 
the in-scattering and the out-scattering contributions to 
the collision integral in the Boltzmann equation which 
emerges from Eq. (A13) in the classical limit. 



2. Distribution function 

Another form of the kinetic equation can be obtained 
by expressing G K in terms of the distribution function g 
by setting 

G K = G R tf -gG A , (A16) 
see Eq. ( 2.26 ). Note that a similar ansatz by Kamenev^H 



assumes that g is hermitian, which in general is not cor- 
rect. A formal proof that G K indeed can be written in 
the form (2.26) can be found in Ref. |48| (see also Ap- 



pendix B of Ref. [17]). The ansatz ( A16| has some formal 
similarity to the generalized Kadanoff-Baym ansatz-^ 

G K {t, t') = i[G R (t, t')G K (t\ t') - G K (t, t)G A (t, t% 

(A17) 

which is often used to express the Keldysh Green func- 
tion at different times in terms of its equal-time matrix 
elements. In fact, the generalized Kadanoff-Baym ansatz 
amounts to setting 



mtt' = [9 f }tt>=6(t-t')iG K (t,t). 



(A18) 



Substituting Eq. ( A16| into our definition (2.25) of the 
Keldysh self-energy E K , we obtain the kinetic equation 



M g - g^M = -t K + t R g - tf± J 



which can also be written as [see Eq. (2.27)] 
- i(Mg - tfM) = E in - E out , 



(A19) 



(A20) 



where E m and E out are defined in Eqs. (2.29 2.30|. Note 



that in the non-interacting limit or in a simple approxi- 
mation where interaction corrections to the K eldvsh self- 
energy Y, K are neglected, we see from Eqs. (2.17a) and 
(2.25) that E K reduces to an infinitesimal regularization 

-ma 



3. Wigner transformed Keldysh Green function 

To describe the approach to equilibrium, it can be ad- 
vantageous to perform a partial Fourier transformation 
of the two-time Green functions with respect to the time 



difference (Wigner transformation), see Eq. (2.33). Using 
the fact that 



[M G K - G K Mojtt 
ro ° dw 
2k 



{id t + id v )G K {t,t') 
e -i"(t-t')id T G K {T;cj), (A21) 



where G K (r; lu) is the Wigner transform of G K (t, t'), we 
obtain the following expression for the Wigner transform 
of our kinetic equation ( A13 ) for the Keldysh Green func- 
tion, 



id T G (t; ui) - [E , G ]( T;w ) - [E , G ]( r;w ) 
= C in {T:Lj)-C out (T;ij), 



where 



G m (r; W ) 
G out (r;a;) 



(r;a;)' 



(t;oi)- 



(A22) 

(A23) 
(A24) 



Recall that the Wigner transform of the product of two 
matrices A and B in the time labels can be expressed in 
terms of the Wigner transforms of the factors A and B 



[AB] {T . u) = / dsi / ds 



X A{t+ —;w 



duii 
ux)B(t 



dco 



2?r 

, S 2 



cj + cj 2 ).(A25) 



Expanding the arguments of A and B for small uj\ , u>2, 
si, S2, and retaining only the first order in the derivatives 
we obtain the approximate expression 

[AB} (T . tU) ^A(r;uj)B(r;u,) 

+ ^ [a T A(T;u)d u B{Tiu)-a u A{T]u)arB(T;u)].{A26) 

Assuming that the r-dependence of the distribution func- 
tions is slow, we may retain only the first term in the 



right-hand side of Eq. (A26), i.e. we approximate the 



Wigner transforms of the products appearing in the ki- 
netic equations by the products of the Wigner transforms 
of the factors. In this approximation the Wigner trans- 
forms of the commu tator terms on the left-hand side of 
the kinetic equation (A22) for G k (t;ui) vanish, while the 



Wigner transforms of the collision integrals are 

C in (r;w) = ^E i^ (T;w)G / (T;u;), (A27) 

C* out (r;w) = iE I {r-uj)G K {T-io). (A28) 

Given G k (t;uj), we still have to perform a frequency 
integration to reconstruct the equal time Keldysh Green 
function (i.e. the distribution function), 

pG K (t;u). (A29) 
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4. Wigner transformed distribution function 

We can avoid the above frequency integration by work- 
ing with the Wigner transform gk(r;u)) of the kinetic 
equation (2.27) for the distribution matrix g, which is 
given in Eq. (2.34). Note that in thermal equilibrium 



and in the absence of interactions the function gk(r;uj) 
is simply given by a symmetrized Bose function 



1 



e /9(w-/i) _ i 



coth 



/3 



(w - fl) 



(A30) 



In an interacting system with well-defined quasi- 
particles, we expect that gk(r;u>) will relax towards a 
t- independent function gk{oj) which for uj close to the 
renormalized quasi-particle energy can be expressed in 
terms of a symmetrized Bose function with some effec- 
tive temperature. Using 

[M g ~ g*M ]tt' = idtQtt' + id t >g*, t - e k {g tv - g* t , t ) 



—iui(t—t') 



2tt 



id T Keg(r; u>) 



-2i(u! — £fc)Img(r; uj) 



(A31) 



we finally arrive at Eq. (2.34) 



APPENDIX B: SECOND ORDER 
SELF-ENERGIES 

For the derivation of the perturbative kinetic equation 



(2.1 ) we need the components of the nonequilibrium self- 



energy to second order in the magnon-phonon coupling. 
In the contour-basis the self-energies are 

K p ' (*. o = «4 £ < W (*> *')Gf_>, o 

1 

+(q^-q)]. (Bl) 

Here the p,p' e {+, — } label the two branches of the 
Keldysh contour. Using the fact that the (q —> —q) term 
gives just a factor of 2 because the Green functions are 
even functions of the momentum argument, we obtain in 
the Keldysh basis 



(t,t>) 



I 

V 



i 

+D^(t,t')G^_ q (t,t')], (B2a) 

q 

+D«(t,t>)Gt q (t,t')], (B2b) 




<--• < •— < + <--• — <- *— < 



G K 



G R 



a c .-■-'^ a<3 a c a' 

< • — < — fc— < + <— #— < — fc— . 

<--« <—♦--< + 



Q 
-< 



9 
< 



D K G A 

Q 



- a Q a 

+ <-4 — <— *— < 

G K 

FIG. 6. (Color online) Diagrams contributing to the magnon 
self-energies to second order in the magnon-phonon coupling. 
The symbols are defined in the caption of Fig. [4] 



The corresponding diagrams are shown in Fig. [6] The 
spectral component of the self-energy is 

q 

+D«(t,t')Gl_ q (t,t')]. (B3) 

Moreover, using the fact that the product of a retarded 
and an advanced function with the same time argument 
vanishes, the Keldysh component of the self-energy can 
alternatively be written as 

s fe (*.*') = ££7^ (*.OG£„(i,f) 
q 

--Dq(*,*') G fe- q (*>*')]• (B4) 

To derive Wigner transformed kinetic equations, we use 
the fact that the Wigner transform of a product of two 
functions with the same time arguments can be written 
as a convolution in frequency space, 



[A(t,t')B(t,t% T ^ 

/OO 
Jse^A{T+-,T--)B{r- 

rli / 

A{t;J)B{t;u-J). 



2tt 



(B5) 



The Wigner transform of the spectral and Keldysh com- 
ponents of our self-energies are therefore given by 



2tT' 



D I J<r,w')GZ_J<r,u>-u') 



L a 



+D^t,t')Gt q (tX)+D^(t,t')Gt q (t,t') . (B2c) 



+D«{T-Lu>)G I k _ q (T-u J -u>)\, (B6) 

O J, / r 

Df( T;w ')G^_,(r;a;-a/) 

,0 Z7T L 

-^(r-.tyjGi.-^u-a/)!. (B7) 
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APPENDIX C: DERIVATION OF THE RATE 
EQUATION 

In this appendix we summarize the approximations 
which are necessary to arrive at the perturbative ki- 
netic equation (2.1 1. Using Eq. (2.26) we may express 



the Keldysh component of the Green function appearing 
on the right-hand side of Eqs. (B6 B7| in terms of the 



Wigner transform g k (j]U)) of the distribution function. 
To leading order in the gradients we may factorize the 
Wigner transform of the matrix products in Eq. (2.26), 



iGk ( r ; w ) ~ G{{t; ui)g k (r; (J). 



(CI) 



After substituting the self-energies (|B6| a nd (|B7[ ) with 
the Keldysh Green function given by Eq. ( |C1[ ) into the 
Wigner transformed kinetic equation (2.38), we obtain an 



integro-differential equation for the nonequilibrium dis- 
tribution function g k (r) = g k (r;u! — e k ). To further 
simplify this equation, let us assume that the phonon sys- 
tem is in thermal equilibrium. Moreover, we approximate 
the phonon propagators by the non-interacting equilib- 
rium propagators. The Wigner transforms of the phonon 
Green functions can then be replaced by the usual Fourier 
transforms, 



where 



' (C2a) 

(C2b) 

D — 0J q — IT] 

i[F*{lo) - F*{u)] = 2nS(u - w q ), (C2c) 



LO — LU q + IT) ' 
1 



1 



1 



Fl( U )=f{w)F£( U ), (C2d) 



coth 



flu 



(C3) 



is the equilibrium phonon distribution. The correspond- 
ing symmetrized phonon propagators are 



1 



(u + irj) 2 - u 2 



D*{u) = \[F*{u) + F%{-u)\ 



(uj — ir)) 2 — uj 2 



(C4a) 



(C4b) 



^h = -[f»-f;h] 

= 7rsgn(w)[(5(a; - oj q ) + S(cu + oj q )], (C4c) 

iDf M = \ [F? (u) + F K q {-u)} (C4d) 

= f{u)D q [uj). (C4e) 



Substituting Eqs. (CI) and (C4e) into Eqs. (B6) and (B7) 



we obtain for the in-scattering and out-scattering com- 
ponents of the nonequilibrium self-energy 



Ejf(r;u;)=i££(T; W ) 



Di(u/)Gi_JT-,u-u/) 



x [f(u')g h - q (T;u-u') + l\ 
b (T;ui) = E£(r;w) fffc (r;w) 



(C5) 



x [/V) + g k - q (r;Lu - J)] g k (T;cj). (C6) 

Combining in- and out scattering contributions to the 
collision integral, we obtain the quantum kinetic equation 

1 /*00 7 / 

x {/V) [9k- q (T;u - u') - g k (T;oj)} 

+1 - 9k-q(r;uJ - u')g k (T;w)\. (C7) 

In thermal equilibrium the different components of 
the magnon Green function satisfy the fluctuation- 
dissipation theorem, which is equivalent with 



3fc (T; W )^coth('^-^) =3°H- 



Using the identity 



coth(x — y) = 



1 — coth x coth y 
coth x — coth y ' 



(C8) 



(C9) 



with x = f3(uj — /i) /2 and y = (3ui' /2 it is easy to see that 
in thermal equilibrium the in- and out-scattering contri- 



butions cancel, so that the right-hand side of Eq. (C7) 
vanishes identically. 



To reduce Eq. (|C7|) to the rate equation (2.1 ) we set 

(CIO) 

~iA . ( C11 ) 



g k {T;u) = 1 + 2n fc (r; w), 
/V) =1 + 2%) = 1 + 



and shift the phonon momentum, q = k — k' . Then 
Eq. ( C7 1 can be written as 



xG k ,(T;u)'){[l + n k (T;uj)]b(uj - d)n k '{T;d) 

-(r;a/)]6(w / -w)]n fc (r;a;)}.(C12) 



n k i 



Finally, we neglect the damping of intermediate states, 
which amounts to approximating 



G£,(t;u/)«27t<Hu/ -e k ,). 



(C13) 
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Then the ^/-integration in Eq. ( |C12 ) is trivial and we 
obtain for 7ife(r) = rifc(r; £&) the kinetic equation 

2 

d T n k (T) = y^2jk-k' D l-k'( e k - <*') 

x |[1 + «fe(r)]6(e fe - e fe ')rz fe /(r) 

-[1 + n fc /(r)][l + 6(e fe - e h /)]n fc (r)}. (C14) 



by Eq. (2.2) 



This is of the form (2.1 ) with transition rates Wk fe' given 
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